Cooling of a lattice granular fluid as an ordering process

We present a new microscopic model of granular medium to study the role of dynamical correlations and the onset of spatial order induced by the inelasticity of the interactions. In spite of its simplicity, it features several different aspects of the rich phenomenology observed in granular materials and allows to make contact with other topics of statistical mechanics such as diffusion processes, domain growth, persistence, aging phenomena. Interestingly, while local observables being controlled by the largest wavelength fluctuations seem to suggest a purely diffusive behavior, the formation of spatially extended structures and topological defects, such as vortices and shocks, reveals a more complex scenario.Comment: 4 pages, 4 figure

Driven granular gases with gravity

We study fluidized granular gases in a stationary state determined by the balance between an external driving and the bulk dissipation. The two considered situations are inspired by recent experiments, where the gravity plays a major role as a driving mechanism: in the first case gravity acts only in one direction and the bottom wall is vibrated, in the second case gravity acts in both directions and no vibrating walls are present. Simulations performed under the molecular chaos assumption show averaged profiles of density, velocity and granular temperature which are in good agreement with the experiments. Moreover we measure the velocity distributions which show strong non-Gaussian behavior, as experiments pointed out, but also density correlations accounting for clustering, at odds with the experimental results. The hydrodynamics of the first model is discussed and an exact solution is found for the density and granular temperature as functions of the distance from the vibrating wall. The limitations of such a solution, in particular in a broad layer near the wall injecting energy, are discussed.Comment: Revised version accepted for publication. New results added and discussions considering tangential forces. 27 pages (19 figures included), to appear in Phys.Rev.

Comment on "Entropy Production and Fluctuation Theorems for Active Matter"

This is a comment to a letter by D. Mandal, K. Klymko and M. R. DeWeese published as Phys. Rev. Lett. 119, 258001 (2017).Comment: 2 pages without figures, in press as a Comment on Physical Review Letter

Self-propulsion against a moving membrane: enhanced accumulation and drag force

Self-propulsion (SP) is a main feature of active particles (AP), such as bacteria or biological micromotors, distinguishing them from passive colloids. A renowned consequence of SP is accumulation at static interfaces, even in the absence of hydrodynamic interactions. Here we address the role of SP in the interaction between AP and a moving semipermeable membrane. In particular, we implement a model of noninteracting AP in a channel crossed by a partially penetrable wall, moving at a constant velocity $c$. With respect to both the cases of passive colloids with $c>0$ and AP with $c=0$, the AP with finite $c$ show enhancement of accumulation in front of the obstacle and experience a largely increased drag force. This effect is understood in terms of an effective potential localised at the interface between particles and membrane, of height proportional to $c\tau/\xi$, where $\tau$ is the AP's re-orientation time and $\xi$ the width characterising the surface's smoothness ($\xi\to 0$ for hard core obstacles). An approximate analytical scheme is able to reproduce the observed density profiles and the measured drag force, in very good agreement with numerical simulations. The effects discussed here can be exploited for automatic selection and filtering of AP with desired parameters.Comment: 13 pages, 3 figure

Growth in systems of vesicles and membranes

We present a theoretical study for the intermediate stages of the growth of membranes and vesicles in supersaturated solutions of amphiphilic molecules. The problem presents important differences with the growth of droplets in the classical theory of Lifshitz-Slyozov-Wagner, because the aggregates are extensive only in two dimensions, but still grow in a three dimensional bath. The balance between curvature and edge energy favours the nucleation of small planar membranes, but as they grow beyond a critical size they close themselves to form vesicles. We obtain a system of coupled equations describing the growth of planar membranes and vesicles, which is solved numerically for different initial conditions. Finally, the range of parameters relevant in experimental situations is discussed.Comment: 13 pages and 5 postscript figures. To appear in Phys. Rev

Driven low density granular mixtures

We study the steady state properties of a 2D granular mixture in the presence of energy driving by employing simple analytical estimates and Direct Simulation Monte Carlo. We adopt two different driving mechanisms: a) a homogeneous heat bath with friction and b) a vibrating boundary (thermal or harmonic) in the presence of gravity. The main findings are: the appearance of two different granular temperatures, one for each species; the existence of overpopulated tails in the velocity distribution functions and of non trivial spatial correlations indicating the spontaneous formation of cluster aggregates. In the case of a fluid subject to gravity and to a vibrating boundary, both densities and temperatures display non uniform profiles along the direction normal to the wall, in particular the temperature profiles are different for the two species while the temperature ratio is almost constant with the height. Finally, we obtained the velocity distributions at different heights and verified the non gaussianity of the resulting distributions.Comment: 19 pages, 12 figures, submitted for publicatio

Noise Rectification and Fluctuations of an Asymmetric Inelastic Piston

We consider a massive inelastic piston, whose opposite faces have different coefficients of restitution, moving under the action of an infinitely dilute gas of hard disks maintained at a fixed temperature. The dynamics of the piston is Markovian and obeys a continuous Master Equation: however, the asymmetry of restitution coefficients induces a violation of detailed balance and a net drift of the piston, as in a Brownian ratchet. Numerical investigations of such non-equilibrium stationary state show that the velocity fluctuations of the piston are symmetric around the mean value only in the limit of large piston mass, while they are strongly asymmetric in the opposite limit. Only taking into account such an asymmetry, i.e. including a third parameter in addition to the mean and the variance of the velocity distribution, it is possible to obtain a satisfactory analytical prediction for the ratchet drift velocity.Comment: 6 pages, 5 figures, to be published on Europhysics Letters; some references have been adde

Non-equilibrium fluctuations in a driven stochastic Lorentz gas

We study the stationary state of a one-dimensional kinetic model where a probe particle is driven by an external field E and collides, elastically or inelastically, with a bath of particles at temperature T. We focus on the stationary distribution of the velocity of the particle, and of two estimates of the total entropy production \Delta s_tot. One is the entropy production of the medium \Delta s_m, which is equal to the energy exchanged with the scatterers, divided by a parameter \theta, coinciding with the particle temperature at E=0. The other is the work W done by the external field, again rescaled by \theta. At small E, a good collapse of the two distributions is found: in this case the two quantities also verify the Fluctuation Relation (FR), indicating that both are good approximations of \Delta s_tot. Differently, for large values of E, the fluctuations of W violate the FR, while \Delta s_m still verifies it.Comment: 6 pages, 4 figure

Critical properties of Ising model on Sierpinski fractals. A finite size scaling analysis approach

The present paper focuses on the order-disorder transition of an Ising model on a self-similar lattice. We present a detailed numerical study, based on the Monte Carlo method in conjunction with the finite size scaling method, of the critical properties of the Ising model on some two dimensional deterministic fractal lattices with different Hausdorff dimensions. Those with finite ramification order do not display ordered phases at any finite temperature, whereas the lattices with infinite connectivity show genuine critical behavior. In particular we considered two Sierpinski carpets constructed using different generators and characterized by Hausdorff dimensions d_H=log 8/log 3 = 1.8927.. and d_H=log 12/log 4 = 1.7924.., respectively. The data show in a clear way the existence of an order-disorder transition at finite temperature in both Sierpinski carpets. By performing several Monte Carlo simulations at different temperatures and on lattices of increasing size in conjunction with a finite size scaling analysis, we were able to determine numerically the critical exponents in each case and to provide an estimate of their errors. Finally we considered the hyperscaling relation and found indications that it holds, if one assumes that the relevant dimension in this case is the Hausdorff dimension of the lattice.Comment: 21 pages, 7 figures; a new section has been added with results for a second fractal; there are other minor change

Which is the temperature of granular systems? A mean field model of free cooling inelastic mixtures

We consider a mean field model describing the free cooling process of a two component granular mixture, a generalization of so called Maxwell model. The cooling is viewed as an ordering process and the scaling behavior is attributed to the presence of an attractive fixed point at $v=0$ for the dynamics. By means of asymptotic analysis of the Boltzmann equation and of numerical simulations we get the following results: 1)we establish the existence of two different partial granular temperatures, one for each component, which violates the Zeroth Law of Thermodynamics; 2) we obtain the scaling form of the two distribution functions; 3) we prove the existence of a continuous spectrum of exponents characterizing the inverse-power law decay of the tails of the velocity, which generalizes the previously reported value 4 for the pure model; 4) we find that the exponents depend on the composition, masses and restitution coefficients of the mixture; 5) we also remark that the reported distributions represent a dynamical realization of those predicted by the Non Extensive Statistical Mechanics, in spite of the fact that ours stem from a purely dynamical approach.Comment: 23 pages, 9 figures. submitted for publicatio